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I.  INTRODUCTION  AND  PERSPECTIVE 


The  published  literature  on  'robust  stabilization  and  performance'  can  be 
viewed  from  two  viewpoints  based  on  the  characterization  of  uncertainty,  namely 
(1)  structured  uncertainty  and  (il)  unstructured  uncertainty.  Much  of  the 
literature  addresses  only  one  of  these  two  types  of  uncertainties  while  there 
exist  few  results  dealing  with  the  combined  uncertainty  case.  The  'structured 
uncertainty'  is  treated  both  in  the  time  domain,  state  space  framework  (where  it 
manifests  itself  in  the  form  of  perturbations  in  the  entries  of  the  matrices 
describing  the  system  behavior)  and  in  the  frequency  domain,  transfer  function 
framework  (where  it  manifests  itself  in  the  form  of  perturbations  in  the  coeffi¬ 
cients  of  the  ncuserator  and  denominator  polynomials  of  the  transfer  function). 
However,  the  'unstructured  uncertainty'  (resulting  mostly  from  the  high  frequency 
unmodeled  dynamics)  is  treated  exclusively  in  the  frequency  domain  framework. 

The  robust  stabilization  problem  in  the  presence  of  structured  uncertainty 
has  been  given  much  attention  in  recent  years.  References  [l]-[23]  provide  a 
representative  list  of  the  type  of  work  being  carried  out  in  this  area.  On  the 
other  hand,  the  robust  stabilization  problem  under  unstructured  uncertainty  has 
been  receiving  attention  much  longer  and  has  produced  many  interesting  results, 
notably  the  H«  theory  and  the  LQG/LTR  theory,  among  others  [24] -[33].  Of  these, 
the  majority  of  the  results  are  concerned  with  stability  robustness  analysis 
whereas  techniques  dealing  with  the  robust  synthesis  problem  are  relatively  few 
[31] -[32].  The  combined  uncertainty  problem,  being  more  difficult,  produced  the 
least  literature  [34] -[39]. 

There  are  three  paths  being  taken  to  address  the  combined  uncertainty 
problem.  One  strategy  is  to  convert  the  structured  uncertainty  into  some  form 
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of  unstructured  uncertainty  formulation  and  then  use  the  methods  available  for 
the  unstructured  uncertainty  problem.  The  structured  singular  value  method  [17] 
falls  into  this  category.  However,  these  methods  are  overly  conservative  because 
they  do  not  take  advantage  of  the  structure  of  the  uncertainty.  The  other  path 
has  been  to  consider  a  'weakly  unstructured  perturbation'  (in  which  we  assume 
both  phase  and  magnitude  information  of  the  uncertainty  are  available,  whereas 
in  'highly  unstructured  uncertainty'  we  assttme  that  only  magnitude  Information 
is  available)  and  then  convert  this  into  a  stnictured  uncertainty  formulation  and 
then  use  the  methods  available  for  structured  uncertainty.  Of  course,  this 
strategy  obviously  is  not  suitable  for  the  combination  of  highly  unstructured 
uncertainty  and  structured  uncertainty.  Finally,  the  third  path  is  to  derive  a 
new  problem  formulation  in  which  the  method  explicitly  accommodates  both  struc¬ 
tured  uncertainty  and  unstructured  uncertainty.  Reference  [37]  was  one  of  the 
first  papers  to  consider  the  design  problem  for  the  combined  case.  In  this 
paper,  this  problem  is  cast  completely  In  the  frequency  domain,  transfer  function 
setting.  Recently,  an  interesting  new  framework  was  proposed  to  solve  the  syn¬ 
thesis  problem  for  the  combined  uncertainty  case  [38] -[39]. 

In  this  research,  some  results  which  consider  only  one  type  of  uncertainty 
as  well  as  some  results  which  treat  the  combined  case  are  presented.  The  report 
is  organized  as  follows:  Section  II  considers  the  robust  stabilization  aspect. 
In  Subsection  2,1,  a  robust  stabilizing  controller  design  for  HIMO  (Multiple 
Input  Multiple  Output)  systems  for  a  class  of  nonminimum  phase  systems  is  pre¬ 
sented.  Subsection  2.2  presents  some  Interesting  results  on  the  simultaneous 
stabilization  of  two  plants  each  having  a  separate  unstructured  uncertainty 
profile.  In  other  words,  this  problem  is  one  form  of  the  combined  uncertainty 
problem  wherein  structured  uncertainty  is  taken  to  be  simultaneous  stabilization 
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(l.e.  the  uncertain  parameter  takes  on  discrete  values  instead  of  continuum  of 
values) .  Subsection  2.3  presents  nonconservative  stability  robustness  bounds  for 
linear  systems  with  real  parameter  variations.  In  Section  III,  the  aspect  of 
performance  robustness  is  addressed.  Finally,  Section  IV  offers  some  concluding 
remarks  and  recommendations  for  future  research. 
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II.  ROBUST  STABILIZATION  IN  THE  PRESENCE  OF  STRUCTURED 
AND/OR  UNSTRUCTURED  UNCERTAINTY 

Here,  we  address  the  issue  of  robust  stabilization  for  linear  systems. 
First,  in  Subsection  2.1,  the  aspect  of  robust  stabilization  under  structured 
uncertainty  (real  parameter  variations)  is  considered.  In  previous  research 
[15],  the  control  design  problem  was  solved  under  some  restrictive  assumptions 
(such  as  minimvun  phase,  etc).  In  this  research,  we  relax  Che  minimum  phase 
assumption,  consider  a  class  of  nonminimum  phase  systems,  and  present  a  control 
design  method  suitable  for  MIMO  systems.  Then  in  Subsection  2.2  the  problem  of 
combined  uncertainty  is  addressed.  Since  it  is  difficult  to  solve  the  control 
design  problem  for  very  general  nonminimum  phase  systems  under  the  combined 
uncertainty  case,  a  different  problem  formulation  in  which  the  structured  uncer¬ 
tainty  is  cast  as  the  simultaneous  stabilization  problem  (for  two  plants)  is 
presented.  In  this  framework,  there  are  no  severe  assumptions  made  on  Che 
nominal  plants.  Finally,  in  Subsection  2.3,  the  problem  of  analyzing  the  sta¬ 
bility  robustness  of  a  given  controller  is  addressed.  This  problem  is  cast  as 
one  of  obtaining  stability  robustness  measures  for  a  given  stable  system.  Using 
known  results  on  the  characteristic  polynomials  of  a  Bialternate  sum  matrix, 
nonconservative  perturbation  bounds  are  obtained  to  maintain  stability. 

2.1  Robust  Stabilization  Under  Structured  Uncertainty 
for  a  Class  of  Nonminimum  Phase  Systems; 

Robust  stabilization  of  linear  uncertain  systems  with  real  parametric 
uncertainty  has  been  an  active  topic  of  research  in  recent  years.  In  fact,  the 
research  documented  in  the  two  final  reports  of  the  contracts  F33615-86-K- 3611 
and  F33615-84-K-3606  carried  out  by  this  principal  investigator  [40] -[41]  serve 
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as  useful  reference  for  this  research  (along  with  references  therein).  In  the 
very  few  methods  available  for  control  design  with  guaranteed  stability  robust¬ 
ness,  a  common  assumption  is  the  minimum  phase  assumption.  Attempts  to  relax 
this  minimum  phase  assumption  and  still  guarantee  robust  stability  turned  out  to 
be  quite  nontrivial,  and  only  partial  inroads  have  been  made  in  this  area. 
Towards  this  direction,  in  [42] ,  a  class  of  nonminimum  phase  Single  Input  Single 
Output  (SISO)  systems  were  considered.  In  this  subsection,  we  present  a  design 
method  useful  for  MIHO  systems,  for  a  class  of  nonminimum  phase  systems,  namely 
those  with  unstable  poles  only  at  the  origin. 

This  problem  is  solved  in  an  indirect  way  as  described  below.  We  first 
discuss  the  robust  stabilization  of  a  ’type  L'  feedback  system.  We  then  show 
that  the  control  design  strategy  for  this  problem  is  directly  applicable  to  the 
robust  stabilization  of  parameter  uncertain  systems  with  unstable  poles  being 
only  at  the  origin  (with  possible  right  half  plane  zeros) . 

2.1.1.  Robust  Stability  of  Tvoe  L  Feedback  Systems: 

In  many  engineering  applications,  control  system  designs  utilizing  some 
type  of  integral  feedback  are  common.  For  example,  it  is  well  known  that  a 
single-input/slngle-output  'type  2'  system  (with  double  integrators  in  the  loop 
transfer  function)  can  track  a  ramp  signal  with  zero  steady- state  error. 
Similarly,  'type  L'  system  design  has  found  use  in  engineering  practice.  In 
[43],  the  concept  of  type  L  system  Is  generalized  to  multivariable  systems. 
However,  it  was  not  until  recently  that  the  robust  stability  of  integral  feedback 
control  ('type  1'  system)  has  received  attention  [44]. 

In  this  subsection,  we  consider  the  robust  stability  of  type  L  feedback 
systems.  We  shall  extend  the  results  reported  in  [44]  to  general  type  L  systems. 
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Unless  stated  otherwise,  we  shall  assume  throughout  the  subsection  chat  the  plant 
Is  an  open  loop  stable,  linear,  time- Invariant  system.  For  simplicity  in 
notation  but  without  loss  of  generality,  we  shall  restrict  the  plant  transfer 
function  to  be  a  square  matrix  relating  p  Inputs  to  p  outputs.  Our  objective  Is 
to  propose  a  simple  design  technique  for  type  L  feedback  systems  which  possess 
certain  stability  robustness  in  the  closed- loop.  Then  this  result  Is  directly 
applied  to  the  problem  of  robust  stabilization  of  parametric  uncertain  systems 
with  right  half  plane  zeros,  a  problem  discussed  In  [42]. 

2.1.2  Preliminaries: 

This  subsection  Is  concerned  with  synthesis  of  type  L  feedback  systems  as 
shown  in  Figure  1.  The  plant  is  assumed  to  be  linear,  time  invariant  and  finite 
dimensional  and  is  described  by  the  following  state-space  model: 


x(t)  -  Ax(t)  +  Bu(t),  y(t)  -  Cx(t) 


(2.1) 


where  x(t)  €  R  ",  u(t)  e  R“,  y(t)  €  R**  with  m  •  p,  and  the  eigenvalues  of  the 
matrix  A  are  all  In  the  open  left  half  plane.  The  transfer  function  matrix  P(s) 
is  represented  by 

P(s)  -  C(sl  -  A)-‘B.  (2.2) 
The  compensator  K(s)  is  assumed  to  be  of  the  form 


Kis)  -  Ki_{a) 


(2.3) 


where  K^'s  e  RP*“,  with  m  -  p  and  L  >  1.  The  prohlra  in  consideration  is  to 
synthesize  the  Kj's  such  that  the  closed- loop  system  admits  certain  stability 
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robustness.  It  wilt  be  shown  that  the  aforementioned  synthesis  task  is  closely 
related  to  singular  perturbation  theory.  The  next  result  will  be  repeatedly 
used  throughout  this  section. 

Theorem  2.1  [45,  p.  52].  Consider  the  singularly  perturbed  (autonomous) 

system 


3t 

>  G 

X 

si 

H  J 

Z 

(2.4) 


where  F  €  R"*",  J  e  and  f^O  is  a  scalar.  Suppose  that  matrices  J  and  F-GJ‘'^H 
are  both  stable  (having  all  eigenvalues  in  the  open  left  half  plane) .  Let 

e*  =  -  ■  (2.5) 

♦  vlSTv'l  ) » 

Then  for  all  «€  (0,  e*) ,  the  system  (2.4)  is  asymptotically  stable. 

Remark  1:  The  estimate  of  e*  in  (2.5)  is  conservative.  The  etact  value 
of  e*  can  be  computed  from  a  more  complex  procedure  as  discussed  In  [46] . 

Using  Theorem  2.1,  we  can  easily  obtain  the  following  result  which  was 
discussed  extensively  in  [44] . 

Theorem  2.2:  Let  the  feedback  compensator  K(s)  be  of  the  form  (2.3)  with 
L  -  1.  Then,  there  exists  Ko  -  toEg  such  that  the  closed-loop  system  is  stable 
whenever  e  €  (0,  Cq)  for  some  €q  >0,  Eq  e  RP*p  and  if  and  only  if  the  matrix 
P(0)  -  -CA'^B  is  nonsingular. 

Proof:  Suppose  the  matrix  CA'^B  is  nonsingular,  then  Eq  e  R““  can  be 

found  such  that  -EgCA'^B  has  all  eigenvalues  in  the  open  left  half  plane.  If  the 
feedback  compensator  as  in  (2.3)  is  used  with  L  -  1,  then  we  have 
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u(t)  -  *0  f  *  fi'oy(t)cft  -  *0  f  *  BoCxlx)dx 
Jo  Jq 


(2.6) 


as  the  feedback  control  law.  Differentiating  both  sides  of  (2.6)  yields 


^d(t)  -  BaCx(t) 


(2.7) 


Ue  may  write  the  above  equation  as  a  unified  matrix  equation  with  (2.1)  as  below: 


O  E, 
B  A 


(2.8) 


Define  «t-r,  then 


dt  dx 


and 


dt  dx 


.  Hence,  (2.8)  can  be  rewritten  as 


du 

dx 


O  EqCI  u(t) 

"  fl  .A  j  U(t)J’ 


(2.9) 


It  should  be  clear  that  the  stability  of  the  closed- loop  system  is  equivalent  to 
the  stability  of  the  singularly  perturbed  system  (2.8).  We  can  now  apply  Theorem 
2.2  to  (2.9)  to  conclude  the  existence  of  xq  >  0  such  that  the  closed- loop  system 
is  stable  for  all  e  €  (0,  « J) .  Conversely,  if  the  closed-loop  system  is  stable 
for  all  e  G  (0,  <o)  for  some  <$  >  0  and  Eg  e  RP***,  then  by  singular  perturbation 
theory,  we  must  have  the  stability  of  matrices  A  and  -EgCA'^B,  i.e.  the  eigen¬ 
values  of  A  and  -EgCA'^B  are  all  in  the  open  left  half  plane.  Note  that  the 
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stability  of  matrix  A  resulting  from  this  theorem  is  superfluous  since  we  have 
already  assumed  it  to  be  stable. 

Clearly,  €q  represents  the  gain  margin  of  the  closed- loop  system.  Hence, 
the  larger  value  of  ej  implies  a  better  stability  robustness. 

Remark  2;  The  value  of  ej  can  be  estimated  from  (2.5)  with  F  -  0,  J  -  A, 
G  -  EqC  and  H  -  B.  Further,  the  matrix  Eq  used  to  maximize  the  estimate 
(2.5).  Further,  Remark  1  can  also  be  used  if  a  maximum  value  of  cq  is  needed. 

It  is  noted  that  the  assumption  on  equal  number  of  inputs  and  outputs  is 
not  a  severe  restriction.  Indeed,  if  the  transfer  function  matrix  of  the  plant 
is  not  square,  then  the  condition  presented  in  Theorem  2.2  can  be  termed  as  the 


matrix 


-A  B\ 

-c  o\ 


having  full  rank.  In  this  case,  we  may  in  principle  employ 
the  technique  in  [47]  to  square  down  the  plant  by  static  or  dynamic  compensation 
while  keeping  the  full  rank  condition.  The  details  are  omitted  here. 


2.1.3.  Main  Result: 

In  this  section,  we  consider  synthesis  of  a  type  L  compensator  as  defined 
in  (2.3)  for  L  >  1.  The  next  theorem  is  the  main  result. 

Theorem  2.3:  Let  the  feedback  compensator  K(s)  be  of  the  form  (2.3)  with 
L  >  1  and  -  «iEi.  Then,  the  closed-loop  system  can  be  stabilized  by  some  K(s) 
of  the  form  (2.3),  if  the  matrix  P(0)  -  -CA'^B  is  nonsingular. 

Proof;  The  theorem  is  true  for  L  -  1  in  light  of  Theorem  2.2.  By 
induction,  it  is  assumed  that  the  theorem  is  also  true  for  L  -  m-1.  That  is,  the 
closed-loop  system  is  stable  for  K(s)  -  Ka,-i(s)  which  is,  in  fact,  equivalent  to 
the  matrix 
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^-1 


A  -BC^, 


(2.10) 


having  all  eigenvalues  In  the  open  left  half  plane  where  (Aj,-!,  6^-1,  Ca-i)  Is  the 

realization  of  K„-i(s)  with  -  «tEi  for  i  -  0.  1 . m-2  and  K„-2  nonsingular. 

We  shall  show  that  >  0  exists  such  that  the  closed- loop  system  is  stable  for 
some  Ea-i  e  whenever  <,-1  €  Indeed,  with  K(s)  defined  as  in  (2.3), 

we  have 


♦  x;.^(s),  (2.11) 


and  the  stability  of  the  closed- loop  system  in  Figure  1  is  equivalent  to  the 
stability  of  the  closed- loop  system  in  Figure  2  because  the  characteristic  equa¬ 
tion  of  the  closed- loop  system  in  Figure  1  can  be  written  as 

det(T  ♦  K.{s)P(3))  ■  detU  ♦  (j  ♦  (s)P(s))‘‘}det(J  +  {s)P(s)  (2.12) 

s" 

E 

Since  K„  i(s)  and  have  the  same  McMillan  degree  and  have  all  the  poles 

TB  i  ffB-X 

at  the  origin,  we  may  always  find  the  realizations  of  K,-i(s)  and  — ^  such 

s 

that  they  have  the  same  A  and  B  matrices.  Let  (A„-i'  ^a-i.  ^m-i)  reali- 

E 

zation  of  Ka.i(s)  and  (A,-!,  B„.i,  C^)  be  the  realization  of  .  Define 

P(a]  m  I^Pia)  {I*K^^{8)P[s))-\  (2.13) 

then  the  realization  of  P(s)  can  be  represented  by  (A,  5,  C,)  with 
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A 

0 

B 

O 

Vi 

0 

,S* 

0 

,  and  d 

0 

0 

^►1 

p 

Cc. 

(2.14) 


Apply  the  similarity  transformation 


T 


I  O  O 
0  10 
O  -I  I 


to  realization  (3.14)  and  eliminate  the  unobservable  states.  We  then  obtain  a 


realization  (A,  B,  C)  of  P(s)  with 


X 


and 


(2.15) 


It  is  noted  that  A  is  the  same  as  Fa-j  as  in  (2.10)  which  is  stable  by  induction. 
Hence,  P(s)  is  stable.  Further,  It  is  noted  that  P(0)  -  Ea-iEa-2  from  (2.13). 
By  Induction,  Ea-2  nonsingular.  This  implies  that  there  exists  Ea-i  e  RP***, 
such  that  all  the  eigenvalues  of  P(0)  are  on  open  left  half  plane.  Therefore, 
if  we  apply  Theorem  2.3  to  the  feedback  system  in  Figure  2  with  P(s)  -  P(s)  and 
cq  “  then  conclude  the  existence  of  >  0,  such  that  the  closed- 
loop  system  is  stable  whenever  €  (0,t*‘^).  Hence,  the  proof  is  complete. 

The  condition  of  nonsingular  P(0)  -  -CA'^B  may  not  be  necessary  for  Theorem 
2.3.  However,  if  we  Intend  to  use  Theorem  2.2  as  a  synthesis  tool  for  the  design 
of  K(s),  then,  nonsingularity  of  P(0)  -  -CA'^B  is  also  necessary.  Next,  we  pre- 
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sent  an  algorithm  for  synthesis  of  K(S)  as  in  (2.3)  based  on  Theorems  2.2  and 
2.3: 

Algorithm  1: 

Let  P(s)  -  C(sl  -  A)"^B  be  defined  as  in  (2.2)  satisfying  the  condition 
that  P(0)  -  -CA'^B  is  nonsingular.  Define  the  type  L  compensator  K(s)  as  in 
(2.3)  with  Ki  -  and  L  i  1.  Then,  can  be  synthesized  as  follows: 

Step  1:  Choose  Eq  €  such  that  -EqCA'^B  has  all  eigenvalues  in  the  open 
left  half  plane; 

Step  2:  Compute  cq  according  to  (2.5)  with  F  -  0,  J  •  A,  G  -  E2C  and 
H  -  B.  Choose  cq  G  (0,  «2)  and  set  Kq  “  «oEo: 

Step  3:  For  m  -  2  to  L,  do  the  following: 

B 

(i)  Define  realizations  for  Kb(s)  and  such  that  they  have 

the  same  A  and  B  matrices.  Denote  the  realization  of  K„(s)  as  (A*,  B,,  ()„)  and 

E 

hence  the  realization  for  — can  be  written  as  (A„,  B^,  Cg)  for  some  C#. 

(ii)  Choose  e  such  that  the  matrix  E„-iEa.2  has  all 

eigenvalues  in  the  open  left  half  plane; 

(iii)  Compute  according  to  (2.5)  with  F  -  0,  J  -  A,  G  -  E^|-iC  and 
H  -  B,  where  (A,  B,  C)  is  defined  as  in  (2.15).  Choose  fn-i  e  (Q.c^*^)  and  set 

Kg-x  “ 

It  should  be  noted  that  in  synthesizing  the  compensator  K(s)  with  Algorithm 
1,  we  may  use  Remark  2  to  achieve  a  better  bound  for  cl's.  However,  it  may 
increase  the  computational  burden. 
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Fig.  2:  Integral  Feedback  Control  System 


13 


Robust  stabilization  of  parametric  uncertain  systems  has  been  studied  exten¬ 
sively  in  recent  years  [48] -[49].  A  common  assumption  used  in  these  references 
is  the  minimum  phase  condition.  Here,  we  intend  to  apply  the  results  in  Subsec¬ 
tions  2.1.2  and  2.1.3  to  a  class  of  parametric  uncertain  systems  without  the 
minimum  phase  condition.  The  uncertain  system  is  assumed  to  have  unstable  poles 
only  at  the  origin  and  have  equal  number  of  inputs  and  outputs  which  can  be 
represented  as 

P(3.q)~-^C(q)  {aI-A{q))-^B(q) ,  (2.16) 

where  Li  1,  and  q  e  Q  is  the  uncertain  parameter  vector.  The  following 
assumptions  are  adopted. 

Assumption  A.l  (continuous  coefficients):  The  elements  of  A(q),  B(q)  and 
C(q)  are  continuous  functions  of  q. 

Assumption  A. 2  (compact  parameter  set):  Q  is  a  compact  set. 

Assumption  A. 3  (stability):  A(q)  has  all  eigenvalues  in  the  open  left  half 
plane  for  all  q  €  Q. 

Assumption  A. 4  (nonsingular  D.C.  gain):  -C(q)A(q)'^B(q)  is  nonsingular  for 
all  q  €  Q  and  there  exists  a  matrix  E  e  such  that  the  eigenvalues  of 
-EC(q)A(q)*^B(q)  are  all  in  the  open  left  half  plane  for  all  q  €  Q. 

Theorem  2.4:  Suppose  the  transfer  function  p(s,  q)  as  in  (2.16)  satisfies 
Assumptions  A1-A4.  Then  there  exists  a  stable  feedback  compensator 

K(s)-N(s)D(s)-^-(Nl-i  +  Nl-jS  +....+  Nos‘'-i)(D,  +  D.-iS  +....+  Dos-'^)'^  (2.17) 

with  m  2  L  such  that  the  closed- loop  is  robustly  stable  for  all  q  €  Q. 

Proof:  Choose  D(s)  “  +  Da.^s  DqS®'^  such  that  det(D(s))  is 
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strictly  Hurwltz  and  det(Do)  »<  0.  Define  P(s,  q)  D(s)‘^C(q)  (si  -  A(q))'^B(q) 
with  realization  (A(q) ,  B(q) ,  C(q)).  Then  for  the  case  L  -  1,  we  may  apply 
Theorem  2.2  to  P(s,  q) .  Indeed,  Assumption  A4  implies  the  existence  of  Eq  e 
such  that  the  eigenvalues  of  -EoC(q)A(q)‘*^B(q)  are  all  in  the  open  left  half 
plane  for  all  q  €  Q.  Hence,  we  may  use  (2.5)  to  estimate  the  value  of  Cg  with 
F  -  0,  J  -  A(q),  G  -  Eo6(q)  and  H  -  6(q).  Since  Q  is  compact  and  A(q) ,  6(q) , 
C(q)  are  continuous  functions  of  q,  «©  >  0  exists.  Hence,  we  may  set  N(s)  -  €oEo 
where  Cq  e  (0,  cq)  .  For  the  case  L  >  1,  we  may  also  follow  the  proof  of  Theorem 
2.3  with  P(s)  replaced  by  P(s,q).  The  same  argument  can  clearly  be  carried 
through  in  the  case  L  >  1  and  hence  the  details  are  omitted. 
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2.2  Simultaneous  Stabilization  Under  Unstructured  Uncertainties: 

As  indicated  in  the  previous  subsection,  the  current  literature  on  the 
robust  stabilization  problem  can  be  attributed  to  treating  two  classes  of  per¬ 
turbations,  namely  i)  parametric  (structured)  uncertainty  and  ii)  high  fre¬ 
quency  unstructured  uncertainty.  In  turn,  the  literature  on  the  structured 
uncertainty  addresses  the  problem  from  two  perspectives:  namely,  the  frequency 
domain  (input/output,  transfer  function)  viewpoint,  and  the  time  domain  (state 
space)  viewpoint.  In  what  follows,  these  developments  are  briefly  reviewed. 

2.2.1.  Parametric  (Structured)  Uncertainty  Problem: 

One  school  of  thought  developed  in  the  area  of  the  robust  stabilization 
problem  is  the  so-called  'simultaneous  stabilization'  concept.  The  objective 
is  to  design  a  controller  which  guarantees  stability  for  a  set  of  plants.  This 
set  of  plants  constitutes  the  uncertain  system  in  which  the  parameters  of  the 
system  are  unknown  but  bounded  within  given  ranges.  The  simultaneous  stabili¬ 
zation  problem  was  introduced  by  Saeks  and  Murray  [50]  and  Vidyasagar  and 
Viswanadham  [51].  Essentially,  the  problem  formulation  is  as  follows:  given  a 
set  of  k  different  plants  {Pi(s),  P2(s),  ...  Pit(s)),  does  there  exist  a  single 
compensator  C(s)  which  stabilizes  the  entire  set?  Clearly  the  problem  of  robust 
stabilization  under  parameter  variations  can  be  cast  in  the  framework  of 
simultaneous  stabilization  as  follows:  Given 

(P  (s,  q)  :  qeQ)  (2.18) 
where  Q  is  a  compact  index  set  and  q  is  an  uncertain  parameter  vector,  does  there 
exist  a  single  compensator  C(s)  to  stabilize  all  the  plants  in  the  set  {P(s,  q))? 
If  the  parameter  vector  q  takes  distinct  (discrete)  values,  then  the  problem  is 
one  of  stabilizing  a  discrete,  finite  number  of  plants  whereas  If  q  takes  on  a 
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continuum  of  values  within  the  bounded  set  Q,  then  It  Is  one  of  stabilizing  an 
Infinite  set  of  plants  (within  the  bounded  set  Q) . 

In  their  paper.  Seeks  and  Murray  [SO]  develop  geometric  conditions  for 
simultaneous  stablllzablllty  but  do  not  offer  any  computational  criterion  for 
Implementing  the  design.  It  Is  mentioned  that  the  computational  criterion  is 
given  only  for  the  two  plant  case  [52].  In  [51]  Vldyasagar  and  Vlswanadham 
generalize  the  above  motion  for  the  multiple  Input  multiple  output  (HIHO)  case. 
Their  result  states  that  the  problem  of  simultaneously  stabilizing  (k  -t-  1)  plants 
Is  equivalent  to  the  problem  of  stabilizing  k  plants  with  the  added  requirement 
that  the  compensator  Itself  be  stable.  In  the  two  plant  case,  this  leads  to  the 
requirement  that  the  difference  plant  (P2(s)  -  P^Cs))  be  stablllzable  via  a 
stable  compensator.  Evidently  one  can  solve  the  2  plant  case  completely  using 
these  results  since.  In  that  case,  the  task  amounts  to  that  of  finding  a  single 
plant  (namely  p2(s)  -  Pi(s))  for  which  one  can  apply  the  results  of  Youla, 
Bonglorno,  and  Lu  [53]  which  essentially  Involves  the  checking  of  the  parity- 
interlacing  property.  In  [51]  It  Is  also  shown  that  ^Iven  two  nxm  plants,  one 
can  genetically  stabilize  them  simultaneously  provided  either  n  or  m  Is  greater 
than  one.  This  result  Is  further  generalized  In  Ghosh  and  Byrnes  [54]  where  It 
Is  shown  that  the  generic  simultaneous  stablllzablllty  of  r  nxm  plants  Is 
guaranteed  If  max  (n,  m)  2  r.  Subsequently  the  problem  of  determining  a 
computationally  feasible  criterion  for  simultaneous  stablllzablllty  Involving 
more  than  2  plants  was  addressed  In  a  paper  by  Emre  [55]  which  presents  a  solu¬ 
tion  for  the  special  case  with  the  Imposition  that  all  the  k  +  1  closed  loop 
systems  possess  the  same  characteristic  equation  which  clearly  limits  the  appli¬ 
cation  of  the  method.  An  Important  result  which  addresses  the  problem  of  simul¬ 
taneous  stabilization  under  a  continuum  of  variations  In  the  plant  parameters  Is 
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the  one  by  Wei  and  Barmish  [15].  In  this  paper,  sufficient  conditions  are  given 
under  which  a  family  of  single  input  single  output  (SISO)  plants  can  be  stabil¬ 
ized  by  a  proper  (or  strictly  proper,  if  desired)  stable  compensator.  Regularity 
conditions  are  imposed  on  the  plant  family  coefficients,  and  it  is  assumed  that 
the  plant  family  is  minimum  phase  with  one  sign  high  frequency  gain.  The  gen¬ 
eralization  of  these  results  to  the  multiple  input  multiple  output  (MIMO)  case 
is  given  in  [16].  In  the  present  research,  these  concepts  are  used  at  the 
problem  formulation  stage. 

2.2.2  Unstructured  Uncertainty  Problem: 

When  the  system  uncertainty  is  mainly  caused  by  high  frequency  unmodeled 
dynamics,  it  can  be  characterized  in  the  form  of  unstructured  uncertainty.  In 
this  framework,  an  uncertain  system  is  represented  by  a  class  of  plants 

M[P„(s),  r(s)]  t  {(I  +  L(s))P„(s):  llL(ja>)ll<|r(Jw)  ]  :V«fco)  (2.19) 
where  Po(s)  is  a  given  nominal  plant,  r(s)  is  a  prespecified  rational  function 
and  L(s)  is  any  unknown  but  stable  rational  matrix  whose  norm  is  bounded  by 
|t(Jct))|  and  the  uncertainty  L(s)  is  a  multiplicative  perturbation.  A  similar 
model  can  be  given  for  an  additive  perturbation  also.  In  this  formulation,  the 
design  model  is  characterized  by  a  given  fixed  nominal  plant  affected  by  a  norm 
bounded  perturbation. 

One  of  the  most  popular  design  methods  in  the  above  formulation  is  the  H* 
method  [24] -[30].  In  this  method,  the  H"  norm  of  the  return  difference  matrix 
of  the  nominal  plant  is  minimized,  thereby  allowing  one  to  find  the  maximum 
tolerable  norm  bounded  (unstructured)  perturbation  (|r(J(j)|)  for  guaranteed 
stability.  On  the  other  hand,  if  the  perturbation  bound  profile  is  given  (i.e. 
if  r(s)  is  known),  then  a  condition  for  maintaining  closed-loop  system  stability 
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is  given  by  Doyle  and  Stein  in  [28]  (see  Chen  and  Desoer  [29]  and  Ridgely  and 
Banda  [33]  for  related  work).  This  condition  essentially  led  to  the  popular 
LQG/LTR  design  technique  [30].  In  a  different  direction,  combining  Nevalinna- 
Pick  theory  and  Youla's  parameterization  of  controllers  [53],  robust  stability 
conditions  are  provided  by  Kimura  [26]  for  single  input  single  output  systems  and 
are  further  generalized  by  Vidyasagar  and  Kimura  [32]  for  multiple  input  multiple 
output  systems.  It  is  this  concept  that  will  be  used  in  this  research  to  address 
the  combined  case. 

2.2.3.  Simultaneous  Stabilization  Under  Unstructured  Uncertainty: 

Consider  the  simultaneous  stabilization  problems  for  the  two>plant  case 
with  each  nominal  plant  having  different  unstructured  uncertainties.  Let  and 
$2  denote  the  sets  of  plants  with  the  nominal  systems  and  Pjo,  respectively. 

Si:  -{Pi:  l|Pi(»-Pio(»!|^|ri(»|}  (2.20a) 
Sj:  -{Pj:  l|P2(»-P2o(»ll.^|f2(»l»  (2.20b) 

where  we  always  assume  P^  and  P2  have  the  same  number  or  fewer  unstable  poles  of 
Pio  and  P20.  respectively.  Suppose  AP:-P2o-Pio  is  strongly  stabilizable ,  i.e.  AP 
satisfies  the  so-called  parity  interlacing  property  (p.i.p.).  The  objective  is 
to  find  a  single  compensator  which  stabilizes  both  S;  and  S2,  simultaneously,  and 
furthermore,  makes  the  closed- loop  systems  satisfy  some  specifications  of  system 
performance . 

At  first,  we  consider  the  relatively  simple  case  where  P^j  itself  is 
stable.  From  the  assumption  commonly  made  for  the  stabilization  under  the 
unstructured  uncertainties,  must  be  stable.  In  this  case,  the  formulation  of 
the  problem  can  be  greatly  simplified.  The  simultaneous  stabilizing  compensators 
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for  Pjo  and  P20  aire  now  parameterized  by  a  free  parameter  R,  say 

C[Pio]  -  (I-RPio)*^R  (2.21) 

provided  R  is  a  stable  stabilizing  compensator  for  the  difference  sys¬ 
tem  dP-P2o-Pio.  When  the  unstructured  uncertainties  are  taken  into  consideration, 
we  at  first  give  a  result  for  robust  stability  due  to  Chen  and  Desoer  [29] . 

Lemma  1;  A  compensator  C  is  a  robust  stabilizer  of  Si  for  the  unstructured 
perturbation  bound  rj,  if  and  only  if: 

1):  C  stabilizes  P^o  and 

li):  I  |C(I+PioC)-^ri|  |.<1 

A  little  modification  will  give  a  condition  for  simultaneous  stabilization 
under  the  unstructured  uncertainties. 

Lemma  2:  a  single  compensator  is  a  simultaneously  robust  stabilizer  for 
and  $2  if  and  only  if  C  stabilizes  Piq  and  P20  simultaneously  and 

I  |C(I+PioC)-iri|  U<1  (2.22a) 

I  |C(I+P2oC)'^r2l  |.<1  (2.22b) 

Now  the  problem  has  been  reduced  to  finding  a  compensator  C€C[Pio]  such 
that  the  norm  constraints  (2.22a)  and  (2.22b)  are  satisfied.  Note  that  if  C  is 
a  simultaneously  stabilizing  compensator  then 

C-(I-RPio)'^R  (2.23) 
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Subsrltute  (2.23)  into  (2.22a)  and  (2.22b),  respectively.  The  robust  stability 
conditions  are  now 


I  |(I-RPio)*^R[I+PPio(I-RPio)'^R]''ri|  U<1  (2.24a) 

I  |(I-RPio)‘^R[I+PP2o(I-RPio)‘^Rr^r2|  |.<1  (2.24b) 

with  R  being  the  strongly  stabilizing  compensator  for  AP.  This  result  can  be 
summarized  as  the  following  theorem. 

Theorem  2.4:  The  two  sets  of  systems  and  S2  can  be  simultaneously 
robustly  stabilized  by  a  compensator  C  if  and  only  if 

||Rri||.<l  (2.25a) 

I  |R(I+APR)-ir2|  U<1  (2.25b) 

where  R  is  a  strong  stabilizer  for  AP.  The  simultaneously  robust  compensator  for 
Sj  and  S2  is  given  by  (2.23). 

The  proof  of  this  theorem  is  available  by  direct  manipulation  of  (2.22a) 
and  (2.22b)  combined  with  (2.23)  and  is,  therefore,  omitted.  From  Kimura  [31] 
we  can  see  our  problem  is  quite  similar  to  the  problem  therein,  say,  an  interpo¬ 
lation  problem  with  strictly  bounded  real  (SBR)  functions  which  Interpolate  the 
given  values  at  some  RHP  points.  Define 

Rti  0 

G(s):  - 

0  R(I>A?R) 
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The  necessary  and  sufficient  conditions  for  simultaneously  robust  stabiliz- 


ability  of  Sj  and  Sj  under  the  unstructured  uncertainties  are 


a)  G(s)  is  a  SBR  function  and 

b)  R  is  a  strongly  stabilizing  compensator  for  AP. 


From  now  on,  we  restrict  our  attention  to  SISO  systems.  Let  us  at  first 
consider  the  coprime  factorization  representation  of  the  systems.  From  [56]  we 
write  the  right  coprime  factorization  (r.c.f.)  representation  of  the  difference 
system  as 

AP  -  ND*^  (2.27) 
where  N  and  D  are  stable  right  coprime  functions.  The  strongly  stabilizing  com¬ 
pensator  for  AP  can  be  now  characterized  by  the  units  in  H*  space  satisfying  some 
interpolation  constraints.  More  precisely,  the  strong  compensator  has  the  form 


R  « 


U-D 

-ir 


(2.28) 


where  U  is  an  arbitrary  unit  which  interpolates  D  at  the  RHP  zeros  of  N.  By 
applying  (2.28)  to  G(s),  the  requirements  for  G(s)  to  be  an  SBR  function  are  now 
equivalent  to 


and 


(2.29a) 


D 

V 


1 


(2.29b) 


Denote  the  zeros  of  N  in  the  RHP  as  Oj,  03..., and  recall  the  definition 
of  the  corresponding  Blaschake  product  B(s) 
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Rz-cV  .  (ai-s)  {aj-s)  .  . .  (aj^-s) 
(a^+s)  (oj+s)  .  .  .  (aj^+s) 


(2.30) 


Since  B(s)  is  all-pass,  (2.29a)  and  (2.29b)  are  now  equivalent  to  the  new 
inequalities  which  read 


(2.31a) 

(2.31b) 

By  denoting 

E,:  - 

(2.32a) 

(2.32b) 

the  interpolation  constraints  are  now  reduced  to  finding  E^eSBR, 

i-1 , 2 ,  which 

interpolate  zero-value  at  a^,  i-1,2 . 1.  Note  that  the  interpolation  problem 

here  is  quite  simple  compared  with  the  standard  Nevanlinna-Fick  Interpolation 
problem  and  possesses  a  very  simple  solution  [57]-[58].  In  fact,  the  desired  SBR 
functions  can  be  parameterized  by  a  free  parameter,  an  arbitrary  SBR  function, 
as 

Ei  -  B(s)Ei,  1-1,2.  (2.33) 

Theorem  2.5:  The  closed- loop  systems  of  S^  and  S2  can  be  robustly  stabi¬ 
lized  by  a  compensator  C  if  and  only  if  two  arbitrary  SBR  functions  E^  and  E2  can 
be  found  such  that 

D*5^€U  (2.34a) 
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(2.34b) 


AP 


where  U  denotes  the  set  of  units  in  H*. 

Proof:  Necessity:  It  is  a  straightforward  manipulation. 

Sufficiency:  Suppose  D+EiNri^— U*  is  a  unit.  From  the  properties  of  the  SBR 
function  and  the  assumptions  for  the  uncertainty  size  r^  and  r2  (without  loss  of 
generality  r^  and  r2  can  always  be  assximed  to  be  a  stable  and  minimal  phase 
function  [6]),  it  is  obvious  that  Eiri^(U*-D)N*^  is  stable.  Take  R*-Eiri^  as  the 
strongly  stabilizing  compensator  for  AP.  By  applying  R*,  (2.25a)  is  satisfied. 
It  can  also  be  verified  that  (2.25b)  is  satisfied  by  substituting  R*  and  (2.34b) 
into  it. 

Since  the  solvability  of  the  problem  in  theorem  2.5  is  not  generally  avail¬ 
able  at  present,  we  can  only  handle  some  simplified  cases  by  making  further 
assumptions  on  the  radius  of  the  unstructured  uncertainties  and  the  difference 
nominal  system  AP.  It  seems  very  possible  to  extend  the  results  to  more  general 
cases  in  future  study. 

Assumption  A5:  The  radii  of  the  unstructured  uncertainties  of  Si  and  $2  are 
constant  and  equal,  l.e. 

ri(s)  -  r2(s)  -  ro 

Assumption  A6:  AP  itself  is  minimal  phase  with  relative  degree  one. 

It  is  clear  with  the  above  assumption  that  any  first  order  unit  which 
interpolates  at  infinity  is  qualified  for  the  stabilization  purpose.  The  first- 
order  interpolation  at  infinity  can  be  simply  reduced  to  some  requirement  for  the 
high  frequency  gain  (HFG)  of  the  desired  unit.  More  specifically,  the  interpola¬ 
tion  constraint  is  satisfied  if  and  only  if  the  HFG  of  the  unit  is  equal  to  the 
Inverse  of  the  HFG  of  AP.  Let  the  HFG  of  AP  be  1/k.  Then  AP  can  be  rewritten 


24 


in  the  form 


1  n  _  1  n  d 
TfS  "  Tc  c  c 


(2.35) 


where  n  and  d  are  monic  polynomials  and  deg(d)-deg(n}-i-l,  c  Is  an  arbitrary  monic 
polynomial  which  serves  as  the  common  determinator  in  forming  the  coprime  factor¬ 
ization  representation  of  AP.  Then  the  desired  unit  should  have  the  form 


''  (s+a, ) 
u  »  ;cn  J ^ 

i.l  (S+t>i) 


(2.36) 


where  a^  and  , i-1 , 2 ,  .  .  . v.  are  all  constants  with  positive  real  parts.  The 
simplest  unit  is  zero-order  which  is  a  constant  k.  It  is  our  observation  that 
in  the  cases  of  minimal  phase  AP,  the  order  of  the  employed  unit  is  not  a  big 
issue.  The  following  is  the  comparison  of  the  cases  with  the  unit  of  order  zero 
and  the  unit  with  order  v.  (2.34)  is  reduced  to 


l«r,[s^]||. 


Ilcr,  I  I.  <  I 


for  unit  of  order  zero  and 


(2.37a) 


(2.37b) 


(c  n  (s*a^)-d  n  is*b^} 


n  n  (s+b,) 
1-1 


.<  1 


(2.38a) 
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krn 


V  V 

c  n  (,s*a.)~d  n  (s+b.) 


n  n  (s+bi ) 
!•! 


dn  (s+b.) 
1-1 


c  n  (s+a.) 
1-1 


for  unit  of  order  V  case.  By  defining 


V 

c:  ■  c  n  (s-t-ai) 
1-1 


•  V 

d;  ■  d  n  (s+bi) 
1-1 


V 

n:  ■  n  n  (s+b.) 
1-1 


we  can  see  the  formulation  of  the  problem  with  (2.38)  is  the  same  as  that  of  the 


problem  with  (2.37).  Now  we  concentrate  our  attention  on  (2.37).  It  will  be 


shown  the  problem  can  be  reduced  to  a  nonlinear  programming  problem  with  linear 


constraints. 


If  c  is  chosen  as  monic  (this  can  always  be  assumed  without  loss  of  gen¬ 


erality)  ,  it  can  be  easily  shown  from  the  definition  of  the  H*  norm  that 


II  d(s) 


€  ^  1 


(2.39) 


From  the  property  of  the  H*  norm  [25],  we  obtain  a  sufficient  condition  for 
simultaneously  robust  stabilization  for  and  $2  under  Assumption  A5  and 
Assumption  A6. 

Theorem  2.6:  and  82-  can  be  simultaneously  robustly  stabilized  by  a 
single  compensator  C  under  Assumption  AS  and  Assumption  A6  if 


ikro 


1 

€ 


(2.40) 
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We  can  choose  the  parameters  of  c(s),  say,  the  zeros  of  c(s),  such  that  the 
H*  norm  of  d(s)c'^(s)  is  equal  to  identity.  Suppose 

m 

d(s)  »  (s^PiXs^pj) .  . .  (s*p,)  -  n  (s*p^) 

ial 

m 

c(s)  »  (s  +  Ci)(s+C2)...(s*^c„)  »  n  (s+c^) 

1-1 

with  Ci,  i-l,2,...,m  are  all  positive.  From  the  multiplicative  property  of  the 
H*  norm,  if  we  choose  Ci>|pi|,  i-1.2 . m  it  always  has 


d(s) 

c(s) 


(s*Pi) 

(S+Cji 


(S*P2) 

(S+Cj) 


(s+ft.)  H 


(S+Pl) 

(S+Ci) 


1 


(2.41) 


Hence,  a  corollary  of  Theorem  2.6  follows  from  the  arguments  between  (2.39)  and 
(2.41). 

Corollary  2.7:  and  $2  can  be  simultaneously  robustly  stabilized  by  a 

single  compensator  C  under  Assumption  AS  and  Assumption  A6  if 


||(cro[£^j||.<l 
under  the  constraints  that 


(2.42a) 


Ci  >  IPil  .  i-1.2, . . .m. 


(2.42b) 


To  verify  the  condition  in  Corollairy  2.7,  it  is  sufficient  to  solve  a 
nonlinear  programming  optimization  problem  under  the  constraints  (2.42b)  with  the 
objective  function  defined  as 


f(c)  -  ||kro 


(2.43) 
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2.3.  Nonconservative  Stability  Robustness  Bounds  Under  Structured 
Uncertainty  for  Linear  State  Space  Models 

The  problem  of  analyzing  the  stability  of  a  family  of  matrices  arises  In 
many  control  systems  applications.  When  the  system  Is  described  by  linear  state 
space  representation,  the  plant  matrix  elements  typically  depend  on  some  uncer¬ 
tain  parameters  which  vary  within  a  given  bounded  interval.  Traditionally,  in 
the  analysis  and  design  of  control  systems,  these  parameters  are  given  a  specific 
nominal  value  and  the  studies  (such  as  stability,  performance,  etc.)  are  carried 
out  using  that  single  value  of  the  parameter  set.  Naturally,  one  Issue  of 
Interest  In  such  a  case  would  be  to  establish  whether  the  system  retains  the 
properties  (such  as  stability)  for  the  entire  given  range  of  parameters  or  not. 
Since  stability  is  usually  the  fundamental  issue,  we  concentrate,  in  this 
research,  on  the  stability  robustness  aspect  of  the  problem.  Consider  the  linear 
state  space  description 

*(t)  -  A  (q)  X  (t)  (2.44) 
where  x(t)€R"  and  q  Is  a  p  vector  of  uncertain  parameters  varying  in  the  pre¬ 
scribed  compact  set  Q.  Specifically,  let  the  parameters  q^  be  given  a  priori 
bounds  as 

ai  ^  qi  ^  qi  (2.45) 
We  can  write  the  matrix  A(q)  of  (2.44)  as 

A(q)  -  Ao  +  E(q)  (2.46) 
where  Aq  Is  an  asjrmptotlcally  stable  nominal  matrix  and  E(q)  is  the  perturbation 
matrix.  This  nominally  asymptotically  stable  matrix  Aq  could  represent  the 
closed  loop  (nominal)  system  matrix  for  a  linear  system.  In  that  case,  the 
problem  of  obtaining  bounds  on  q^  to  maintain  stability  is  indirectly  related  to 
the  assessment  of  stability  robustness  of  the  control  gain  matrix. 

The  above  issue  of  determining  the  stability  of  an  Interval  parameter  matrix 
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has  attracted  considerable  amount  of  research  in  the  last  few  years.  A.  very 
Informative  account  of  this  research  is  elegantly  summarized  in  refs.  [11-[21 
(and  their  bibliographies) . 

There  are  essentially  three  paths  being  followed  to  answer  the  question  of 
interval  pararaete'^  matrix  stability. 

One  framework  is  based  on  the  polynomial  theory,  where  the  matrix  stability 
problem  is  converted  to  a  polynomial  root  testing  problem  and  the  results  avail¬ 
able  in  the  uncertain  polynomial  theory  are  used.  The  references  [3] -[4]  belong 
to  this  category.  The  second  viewpoint  is  the  direct  matrix  theory.  That  is, 
trying  to  solve  the  problem  without  converting  to  polynomials  using  such  time 
domain  techniques  as  Lyapunov  theory  and  other  matrix  domain  results.  The 
references  [5] -[14]  fall  into  this  category.  Finally  the  third  route  follows  the 
frequency  domain  (transfer  function)  viewpoint  wherein  the  concepts  of  multi - 
variable  stability  margin  (MSM)  and  structured  singular  value  (SSV)  are  employed 
(as  measures  of  stability  robustness).  The  research  of  references  [17] -[22] 
adopts  this  route.  Host  of  the  results  to  date  on  this  problem  are  essentially 
in  the  form  of  sufficient  conditions  except  for  these  references.  The  results 
of  [3]  consider  special  cases  such  as  n-2  with  only  real  spectrum.  The  recent 
paper  of  [20]  uses  an  iterative  perturbation  domain  splitting  technique  to  close 
the  gap  between  the  necessary  bound  and  the  sufficient  bound.  In  this  research, 
we  present  highly  nonconservative  bounds  on  the  parameters  to  maintain  stability. 
For  some  problems,  these  are  almost  necessary  and  sufficient  bounds. 

In  the  next  subsection,  we  recall  a  result  that  is  fundamental  to  the  devel¬ 
opment  of  the  robustness  problem,  and  in  Section  III  we  present  the  main  result. 
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In  Jury  [59] ,  the  following  necessary  and  sufficient  condition  Is  stated  for 
the  stability  of  an  nxn  real  matrix  A<}  in  terms  of  the  elements  of  the  matrix  Ag 
which  was  originally  due  to  Fuller  [60].  We  reproduce  the  theorem  as  It  Is 
stated  in  Jury  [59]. 

Theorem  2 . 8  [59]:  Let  -  [a^j]  be  a  real  square  matrix  of  dimension  n  >  1. 
Let  Ao-[apq  r»]  'J®  the  'blaltemate  sum'  of  Ao  Itself.  That  Is,  let 

Ao  -  2Ao*I„  -  A*I„+VA  (2.47) 
where  •  denotes  blaltemate  product.  Thus,  A^  Is  a  square  matrix  of  dimension 
m  -  hn(n-l)  with  rows  pq  (p-2,3,...n;  q-1,2. . .p-1) ,  columns  rs  (r-2,3,.,.n, 
s-l,2...r-l)  and  elements  given  by 

**p« 


If  r-q 

if  r-p  and  s-q 

if  r-p  and  s-q  (2.48) 

if  r-p  and  s-q 

if  s-p 
otherwise 


For  the  characteristic  roots  of  Ao  to  have  all  their  real  parts  negative  (i.e. 
stable  matrix),  It  Is  necessary  and  sufficient  that  In  (-1)"  times  the  char¬ 
acteristic  polynomial  of  Ag,  namely 

(-1)"  |Ao-AI„|  (2.49) 

and  In  (-1)"  times  the  characteristic  polynomial  of  Ag,  namely 

(-!)■  |Ao-mI,|  (2.50) 

the  coefficients  of  (1-0, 1,  .  .  .n-1)  and  -  (A^+AJ)^  (1-0,  l,...m-l)  should  all 
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be  positive. 

Before  going  into  the  implications  of  this  theorem,  we  borrow  some  examples 
from  [59]  to  illustrate  the  form  of  for  a  few  matrices. 

For  n-2 


For  n-3 


^11  ®12 
®21  *22 


^  ■  "(*11  *22) 


(2.51a) 


*11  *22  *13 
*21  *22  *33 
*31  *32  *33 


For  n-4 


^0 


ll'*'*22 

*23 

"*13 

*32 

*33**11 

*12 

-*31 

*21 

*33**2 

*11 

*12 

*13 

*14 

*21 

*22 

*23 

*24 

*31 

*32 

*33 

*34 

**l 

*42 

*43 

*44 

>11**22 

*23 

"*13 

*24 

"*14 

0 

*32 

*11**33 

*12 

*34 

0 

"*14 

"*31 

*21 

*22**33 

0 

*34 

"*24 

*42 

*43 

0 

*11**44 

*12 

*13 

-*41 

0 

*43 

*21 

*22**44 

*23 

0 

"*41 

"*42 

*31 

*32 

*33**44 

(2.51b) 


(2.51c) 


In  [59]  ,  a  method  to  form  Aq  from  A^  is  presented  and  for  brevity  reasons, 
this  method  is  not  discussed  here.  It  suffices  to  mention  that  it  is  computer 
amenable  and  easy  to  build. 

Evidently,  the  above  theorem  is  a  very  powerful  and  useful  theorem  on  the 
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stability  of  the  matrix  Ao.  The  main  reason  that  this  theorem  forms  the  backbone 
of  our  further  analysis  on  uncertain  matrices  is  the  property  of  'convexity'  in 
the  individual  coefficients  (of  the  characteristic  polynomials  of  Aq  &  Aq)  to 
guarantee  stability  of  Aq.  By  constructing  the  matrix  Ao,  which  henceforth  will 
be  called  the  'tilde*  matrix  instead  of  the  long  name  'blalternate  sum  of  Aq  with 
itself'  ,  we  can  do  away  with  the  interdependence  of  the  coefficients  (in  the  form 
of  Hurwitz  conditions).  Of  course,  the  price  we  pay  for  this  simple  'positivity 
testing'  condition  is  the  construction  of  Ao  which  is  of  higher  dimension. 
Recall  the  well  known  fact  that  the  Hurwitz  region  of  2x2  matrices  is  convex  in 
coefficient  space.  The  above  theorem  is,  in  a  way,  a  generalization  of  this  fact 
for  nxn  matrices.  Clearly  this  important  property  of  convexity  in  individual 
coefficients  proves  to  be  extremely  useful  in  our  extension  to  uncertain 
matrices . 

It  may  be  recalled  from  (2.45)  that  the  perturbation  parameter  domain  is  a 
p  dimensional  'hyperrectangle'  with  r-2*’ vertices  whose  outer  (exposed)  edges  are 
parallel  to  the  coordinate  axes  q^.  However,  in  sequel,  it  proves  to  be  useful 
to  work  with  a  symmetric  range  in  the  parameters  and  hence  we  will  describe  the 
interval  parameter  matrix  (2.44)  in  an  alternate  form  as  follows: 

The  interval  parameter  matrix 

A  (a)  -  Ao  + 

can  be  written  as 

M  (£)  -  M*  +  T  ejAi  (2.52) 

i-i 

where  the  uncertain  parameters  e^  have  a  symmetrical  range  of  [-Gi,  Gj  ,  tGi>0]  , 
about  the  'center'  matrix  M*  -  M(0) .  Note  that  Gj,  and  M*  can  be  determined 
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easily  from  the  bounds  and  the  vertex  matrices  M|. 

Let  €j.  Then  k^j  determines  the  shape  of  the  hyperrec tangle  In  the 

parameter  space. 

We  now  state  an  obvious  theorem  on  Interval  parameter  matrix  stability  as 
follows : 

Theorem  2.9:  The  Interval  parameter  matrix  (2.52)  Is  stable  If  and  only  If,  In 
(-1)"  |M  (£)  -  A  I„|  (2.53a) 

and  (1)“  I  ft  (£)  -  /i  I„|  (2.53b) 

the  coefficients  of  (which  we  will  call  1-0,1,  ...  n-1)  and  the  coeffi¬ 
cients  of  (which  we  call  8^,  1-0,1,  ...  m-1)  are  all  positive. 

Let  -  {^1  Mi)  (2.54) 

However,  this  theorem  Is  hardly  useful  from  a  practical  viewpoint.  What  we 
need  Is  a  finitely  computable  test  or  condition  to  establish  the  stability  of  all 
the  matrices  In  the  hyperrectangle.  With  this  In  mind,  we  present  a  sufficient 
test  for  establishing  the  stability  of  the  uncertain  matrix  (2.52). 

Theorem  2.10:  The  interval  parameter  matrix  of  (2.52)  is  stable  if 

Max  €i  -  e  <  Mn  (2.55) 

where  Mn  obtained  as  follows: 

1)  Obtain  the  coefficients  7i  of  (2.54)  In  terms  of  the  parameters  e^ 
11)  Obtain  the  polynomials  7i  (of  ei)  where  the  coefficients  of  the  terms 
In  the  polynomials  are  all  made  negative  (except  for  the  constant  term 
which  Is  always  positive  since  It  corresponds  to  the  stable  matrix 
M*). 
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lil)  Obtain  the  single  variable  polynomials  7^  from  the  multivariate 
polynomials  7^  by  substituting  -  k^j  Ej 

Iv)  Solve  for  the  positive  real  roots  of  all  these  single  variable 
polynomials . 

v)  Then  -  minimum  of  all  the  positive  real  roots. 

Proof  :  For  the  interval  parameter  matrix  M(e)  of  (2.52),  the  coefficient  7^  of 
the  characteristic  polynomials  of  M(e)  and  H(e)  are  polynomic  in  the  variable  e^. 
For  example,  when  there  are  two  uncertain  parameters  ej  and  62,  a  typical 
expression  for  the  coefficient  7^  (say  73)  is  of  the  form: 

73  ”  b(j3  +  bj^3e^  4-  b23e2  +  b33e^e2  +  b43ef  +  b33e2  +  ^530^ 

+  ^2302  +  bg3eje2  +  bg3eie2 

where  the  constant  bo3  is  positive  (corresponding  to  the  stable  matrix  M*) 

Now  letting  jeil  -  e^^  and  |e2|  -  ^e  write 

73  -  bo3  -  |bi3|eia  -  |b23|e2,  -  |b33|ei„e2„  -  [b^jle?* 

l^53l®2m  ■  |b83l®L  *  l**73l®L  *  1^831®!*®!® 

-  1^93  I  eL  ®2« 

Then  using  e^^  -  ki2  ®^  where  ki2  >  0  is  known,  we  get 

73  ”  bo3  -  |bi3|ki2®2B  -  1^23 1  e2„  -  1^331^12  ®aB 

■  |b♦3|kf2  ®iB  -  |bj3|e|,  -  |bg3|ki^  e|, 

-  |b73|e2,  -  |bg3|ki2  -  jbgslkf*  e^ 

Clearly,  this  is  a  single  variable  polynomial  in  e2i„. 

We,  then,  solve  for  the  positive  real  roots  of  73  -  0  and  denote  the  minimum 
of  these  as  a3.  If  there  are  no  positive  real  roots  of  this  equation,  it  means 
«3  -  Clearly  73  >  0  for  all  ej*  <  <*3  and  73  >  0  for  all  e2a  -  |e2|  <  <*3  and  eja, 

-  |ei|  <  ki2  «3. 
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Clearly  -  Min 

Note  that  the  proposed  sufficient  condition  is  valid  for  any  general  type 
of  variation  of  the  parameter  vector  in  the  matrix  E(q) . 
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III.  PERFORMANCE  ROBUSTNESS  IN  LINEAR  UNCERTAIN  SYSTEMS 


As  mentioned  earlier,  much  of  the  published  literature  on  robustness 
essentially  deals  with  the  stability  robustness  aspect.  However,  it  is  very  well 
recognized  that  present  day  control  systems  are  required  to  not  only  stabilize 
the  plant  but  also  achieve  some  prescribed  level  of  performance.  Thus,  perform¬ 
ance  robustness  is  an  Important  feature  of  any  controller  design  for  realistic 
application  problems.  Studies  discussing  the  performance  robustness  aspect  have 
been  relatively  scarce  since  it  is  a  much  more  difficult  problem.  Here  we  use 
the  term  'performance'  in  a  qualitative  manner,  which  may  include  achieving  low 
regulation  cost  with  low  control  effort  in  regulation  problems  or  achieving 
desired  transient  and  steady-state  dynamic  response  with  minimias  percentage  over¬ 
shoot,  damping  ratios  and  steady- state  errors,  etc.  Yet  another  way  of  measuring 
performance  is  disturbance  rejection.  It  is  well  recognized  in  eigenstructure 
assignment  theory  that  the  majority  of  these  performance  constraints  can  be 
dequately  modeled  in  terms  of  desired  closed  loop  pole  regions  in  the  complex 
plane.  Hence,  in  this  research,  the  performance  robustness  problem  is  cast  as 
a  ' D-stablllty'  problem  where  the  'D-regions'  are  regions  in  the  complex  plane 
(which  are  symmetric  with  respect  to  the  real  axis  and  are  simply  connected)  in 
which  the  desired  closed  loop  poles  are  assigned.  Thus,  assuming  the  'nominal' 
system  achieves  O-stablllty,  the  performance  robustness  problem  is  taken  to  be 
the  problem  of  achieving  D-stability  under  perturbations.  Another  reason  the  D- 
stability  problem  formulation  is  preferred  is  that  it  becomes  amenable  to  treat 
discrete  system  stability  (where  the  D- region  is  a  unit  circle  in  the  complex 
plane) . 

In  the  proposed  D-stability  problem  formulation  the  uncertainty  under 
consideration  is  termed  'weakly  unstructured.'  As  mentioned  earlier,  'highly 
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unstructured  uncertainty*  is  meant  to  be  the  case  where  only  a  magnitude  bound 
on  the  norm  of  the  perturbation  is  assumed  known.  In  this  case,  the  magnitude 
bounding  function  may  not  be  realized  by  a  rational  transfer  function.  On  the 
other  hand,  in  the  weakly  unstructured  case,  the  bounding  function  can  be 
realized  by  a  finite  dimensional  rational  transfer  function.  The  0-stability 
problem  under  a  highly  unstructured  uncertainty  problem  is  not  amenable  for 
analytical  formulation.  Hence,  in  this  research,  the  D- stability  problem  is 
formulated  for  weakly  unstructured  uncertainty. 

3.1  D-Stablllty  Problem  Under  Weakly  Unstructured  Uncertainty 

In  what  follows,  we  make  these  assumptions  on  the  D-region: 

The  D-region  is: 

1):  Symmetric  with  respect  to  the  real  axis, 

il):  Simply  connected 
Lemma  3.1:  [SISO  Systems]: 

Suppose  g(s)  is  a  rational  continuous  transfer  function.  If  g(s)  has  7 
poles  outside  DUaO  [i.e.  i-1...7.  D  is  complementary  of  D]  then 

g(s)/(l+g(s)  has  p+7  poles  in  D  if  "-1"  Is  not  on  the  Nyquist  locus  along  30,  and 
it  encircles  "-1"  p  times  in  clockwise  sense. 

Lemma  3.2: 

The  closed- loop  system  g(s)/(l+g(s)  is  0-stable  if  the  Nyquist  locus  along 
30  encircles  "-1"  -7  times  in  the  clockwise  sense. 

Remark:  D-stabllity  of  g(s)/(l+g(s)  la  assured  if  and  only  if 

|l+g<s)|,.B  "  0  or  |l+g(s)|„5>0  (3.1) 

Lemma  3.3:  [Robust  0-stabillty] : 

Suppose:  i)  g(s)  -  go(s)  +  dg(s) 
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ii)  go(s)k(s)/(l+go{s)k(s),)  Is  D-stable. 

Ill)  Ag(s)  Is  D-stable.  and  ||Ag»R‘^||.  <  1 
Then  the  closed-loop  system  g(s)k(s)/(l+g(s)k(s) )  Is  D-stable  If: 

|l+(go(s)  +  Ag(s))k(s)  |,,,D  0  (or.  >0) 

Proof:  Since  the  Nyqulst  locus  r[l+go(s)k(s) ]  encircles  the  origin  7  times, 

and  the  number  of  D-unstable  poles  of  go(s)k(s)  and  of  [go(s)+Ag(s)  ]k(s)  are  the 
same  (l.e.  Ag(s)  Is  D-stable),  the  closed-loop  system  g(s)k(s)/(l+g(s)k(s) )  Is 
robust  D-stable  If  r(l+(go(s)+Ag(s))k(s))  encircles  the  origin  -7  times.  Now, 
we  need  to  show  that: 

Number  of  encirclements  of  r(l+(go+Ag)k)  about  "0" 

-  Number  of  encirclements  of  r(l+gok)  about  "0"  * 

-  -7 

This  Is  true  If  and  only  If 

|l+(8o(s)  +  <*Ag(s))k(s) ^  0  (or  positive)  for  all  0  S  e  S  1 
Sufficiency: 

Suppose  |l+(go(s*)  +  <«Ag(s*))|  -  0  where  s*£9D,  the  number  of 
encirclements  of  r(l+(g(,(s)  +  fAg(s)k(s))  about  "0"  r*  -7.  Since  r(*) 
Is  continuous  with  e,  this  contradicts  with  statement  (*) . 
Necessity: 

Since  l+(gQ+<«Ag)k  Is  a  continuous  function  of  « ,  If  ()  Is  not  true, 
there  must  be  a  05e*il  such  that,  r(l+(go+c*Ag)k)  passes  through  "0"  . 
l.e.  I l+(go(s*)+«*Ag(s*) I  -  0,  s*e3D.  Therefore,  a  contradiction. 
This  completes  the  proof. 

Remark:  We  have  the  following  equivalence: 

|l+(go+<*Ag)k|„ao  >  0 

If  and  only  If 
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I  (1+gok)  (l+€Ag»k(l+gok)'M  l.«*D  >  0 
If  and  only  If 
(l+«Agk(l+gok)-»|  ■caD  ^  C 
if  and  only  if 

|l+Agk(l+gok)-i| 

■(to  ^  ® 

Nov  we  are  in  a  position  to  state  the  following  theorem: 

Theorem:  3.1  [Robust  D-stability] : 

The  closed-loop  system  g(s)k(s)/(l-(-g(s)k(s)) ,  where  g(s)-go(s)+Ag(s) ,  is 
D- stable  if: 

|Ag(s)k(s)(l+go(s)k(s)'^)  [■.ao  <  1  (3.2) 

In  order  to  use  the  above  condition  for  design  purposes,  we  take  note  of  the 
following  result: 

Lemma  3.4:  [61] : 

Every  simply  connected  region  D  in  the  complex  plane  is  conformally  equiva¬ 
lent  to  the  open  unit  disc  U. 

Remark:  Obviously,  D  is  also  conformally  equivalent  to  the  LHP. 

Lemma  3.5: 

Suppose  4  is  the  conformal  mapping  function  D-»LHP.  For  any  real 

function  g(s) ,  we  have: 

Sup  |g(s)(  -  Sup  |g(^*^(z))| 
s-jw  ZeflD 

or:  Sup  |g(z)|  -  Sup  |g[^(s)]| 

ZedD  s-Jw 

Define  Hq  as  the  following  Banach  space: 

Ho  t  (f(s):  f(s)  is  analytic  in  D  and  Sup.  |f(s)  is  bounded). 

scD 

and  I |f(s)| lo  i  Sup  |f(s)| 

SeD 
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Remark:  If  D  LHP  and  g<s)  is  D-scable  then  ||g(s)||,  5  ||g(s)||o.  This  is 

called  the  Maximum  Modulus  Theorem  [61], 

The  robust  D- stability  condition  is  given  by 

||Agk(l+gok)-MlD  <  1 
or: 

I  |Ag(^(s))k(^(s))(l+go(^(s).k{^(s))]-M  I.  <  1 
Denote  f (• )  -  k(. ) [l+go(‘ )k(. )  1'^ 

We  reduce  the  condition  to: 

||Ag(d(s)).f(^(s))|U  <  1  (3.3) 

However,  in  practice  we  only  know  the  frequency  dependent  bounding  function  on 
the  uncertainty,  i.e.,  |Ag(jw)|  <  jRCjw)).  The  question  then  is  to  be  able  to 
get  a  condition  for  D> stability,  knowing  only  R(Jw)  (i.e.  ,  frequencies  only  along 
the  Imaginary  axis) . 

ThS!,9g.SB.J .  2 : 

Suppose  we  can  make  an  estimation 

|Ago^(jw)/Ag(Jw)|  <  |tf(jw)|  (3.^) 

Then:  the  closed- loop  system  is  D- stable  if: 

lf(^(Jw))-R(Jw).S(jw)|  <  I 
Proof:  (3.4)  can  be  rewritten  as: 

I  4g(^(Jw))»f(^(Jw))«Ag(Jw)j;~Hjw)  ,  .3  gv 

Ag(Jiir)«R*^Jw) 


(3.5)  can  be  guaranteed  if:  ( |Ag<R'^| )<1. 

i.e.  |f(^(jw)).R(Jw)^(Jw)|  <  1  (3.6) 
If  Ag  is  strictly  unstructured,  the  estimation  of  | Ag(^(Jw) )/Ag(Jw) |  may  not 
exist.  Consider  the  weakly  unstructured  case,  i.e.  when  Ag(s)  can  be  written  as 
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a  rational  function  with  paraneter  variation: 


Ag(8)  -  E(s,q)  qeQ 

In  this  case,  an  estimation  of  |  ^  obtained  and  finally  the 

robust  D-stability  problem  can  be  reduced  to  an  ordinary  robust  stability  problem 
*  which  may  be  solved  using  the  available  techniques. 

Consider  a  simple  case: 


Straightforward  calculation  gives: 

ls*(w)5  -1.618  (3.8) 

.  6^/r -18 

In  this  case  tf(jw)  -  1.618. 

A  more  general  case  is  given  as  follows , 
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Direct  calculation  gives: 


fmax  J 


JcV4jc^»l  -4k*-k^ 

,  (3k*+l)A^c^+l  -4k*-5k*-l 


(3  10) 


The  plot  fmax  (k)~k  is  given  In  fig  3  to  show  the  significance  of  k  in  estimation 
of  4(w) .  It  can  be  seen  when  k  >  0.2 

fmax  <  5.8. 

and:  k  -►  *,  fmax  -»  1. 

Remark:  There  Is  one  important  point  to  be  mentioned  here;  the  final  estimation 
of  4(w)  is  not  related  to  the  uncertain  parameter  bounds .  Of  course  this  is  only 
true  for  first  order  transfer  functions. 
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VI.  CONCLUSIONS  AND  RECOMMENDATIONS  FOR  FUTURE  RESEARCH 
4.1  Work  In  Retrospect: 

The  main  theme  of  the  described  research  under  the  present  contract  has 
been  to  analyze  and  synthesize  controllers  for  robust  stability  and  performance 
under  the  presence  of  structured  and  unstructured  uncertainties.  First  the 
aspect  of  robust  stabilization  under  structured  uncertainty  is  considered.  A 
class  of  nonmlnlmvuo  phase  MIMO  systems  with  structured  uncertainty  is  considered 
and  a  robust  control  design  scheme  Is  presented  using  the  concept  of  Integral 
feedback  control.  Then  the  aspect  of  simultaneous  stabilization  (of  two  plants) 
under  unstructured  uncertainty  is  considered,  and  a  method  for  determining  a  sta¬ 
bilizing  controller  using  concepts  is  presented.  These  two  sections  cover  the 
task  described  as  'task  1'  in  the  original  proposal. 

Next,  the  aspect  of  analyzing  the  stability  robustness  of  linear  systems 
under  structured  uncertainty  is  addressed.  In  this  research,  the  state  space 
description  of  the  linear  system  subject  to  real  parameter  variations  is  taken 
as  the  uncertain  system.  Using  the  theory  of  characteristic  polynomials  of  some 
specially  constructed  matrices,  a  method  for  obtaining  nonconservative  stability 
robustness  measures  is  presented.  This  section  covers  the  task  defined  as  'task 
2'  of  the  original  research. 

Finally  the  aspect  of  performance  robustness  in  linear  uncertain  systems 
is  addressed.  The  performance  robustness  problem  is  cast  as  a  'D-stabilization' 
problem,  and  a  method  that  leads  to  the  design  of  controllers  to  keep  the  closed- 
loop  poles  in  a  designed  region  under  'weakly  unstructured'  uncertainty  is 
presented.  This  result  covers  the  term  labeled  as  'task  3'  of  the  original 
proposal . 

As  it  normally  occurs  in  an  open  research  effort  such  as  this  one,  the  time 
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schedule  for  carrying  out  these  tasks  did  not  adequately  match  the  time  schedule 
originally  proposed  in  the  proposal.  The  research  carried  out  under  task  2  con¬ 
sumed  much  more  time  than  anticipated. 

The  publications  listed  as  references  [62-63-64]  are  the  result  of  this 

study. 

Topics  for  Further  Research 

>■  1)  An  interesting  area  of  research  is  in  the  lines  of  the  work  being  carried 

out  in  [39]  ,  labeled  'mixed  H2/H«'  problem.  This  problem  formulation  seems 
amenable  to  carry  out  robust  stability  and  performance  studies  under  the 
combined  (simultaneous  presence)  uncertainty  case. 

2)  Another  topic  for  further  research  would  be  to  obtain  necessary  and  suf¬ 
ficient  conditions  for  robust  stability  of  linear  interval  parameter 
systems . 

3)  One  research  area  that  needs  attention  is  simultaneous  ei gens true ture 
(eigenvalues  as  well  as  eigenvectors)  assignment  for  linear  uncertain 


systems . 
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